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The observed suppression of high-pT hadron spectra, finite azimuthal anisotropy, disappearance of 
jet-like back-to-back correlations, and their centrality dependence in Au + Au collisions at RHIC are 
shown to be quantitatively described by jet quenching within a pQCD parton model. The difference 
between h and n° suppression in intermediate pr is consistent with the observed (K + p)/ir 
enhancement which should disappear at pr > 6 GeV/c. The suppression of back-to-back correlations 
is shown to be directly related to the medium modification of jet fragmentation functions (FF) similar 
to direct-photon triggered FF's. 
PACS numbers: 12.38. Mh, 24.85. +p; 25.75.-q 
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The degradation of high-p^ partons during their prop- 
agation in the dense medium can provide critical infor- 
mation necessary for detection and characterization of 
the strongly interacting matter produced in high-energy 
heavy-ion collisions. Because of radiative parton en- 
ergy loss induced by multiple scattering, the final high- 
Pt hadron spectra from jet fragmentation are expected 
to be significantly suppressed [1]. Such a phenomenon, 
known as jet quenching, was observed for the first time in 
Au + Au collisions at the Relativistic Heavy-ion Collider 
(RHIC) [2,3]. One also observes the disappearance of 
back-to-back jet-like hadron correlations [4] and finite az- 
imuthal anisotropy [5] of high-pj 1 hadron spectra. These 
three seemingly unrelated high-p^ phenomena are all pre- 
dicted as consequences of jet quenching [1,6-8]. Together, 
they can provide unprecedented information on the prop- 
erties of dense matter produced at RHIC. 

In this Letter, we will study these three high-p^ 
phenomena simultaneously within a lowest order (LO) 
pQCD parton model that includes initial nuclear Jct 
broadening, parton shadowing and medium induced par- 
ton energy loss. We point out that an enhanced (K+p) /ir 
ratio leads naturally to different suppression of /i ± and n° 
spectra at intermediate pt range. We will also show that 
the suppression of back-to-back correlations is directly 
related to the medium modification of hadron-triggered 
FF's similar to a direct-photon triggered FF [10]. 

In a LO pQCD model [9] , the inclusive high-p-p hadron 
cross section in A + A collisions is given by 
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tering cross sections and £a(&) is the nuclear thickness 
function normalized to J d 2 btA(b) — A. We will use a 
hard-sphere model of nuclear distribution in this paper. 
The K w 1.3 — 2 factor is used to account for higher order 
pQCD corrections. The parton distributions per nucleon 
fa/A( x aiQ 2 T r ) inside the nucleus are assumed to be fac- 
torizable into the parton distributions in a free nucleon 
given by the MRS D— ' parameterization [11] and the 
impact-parameter dependent nuclear modification factor 
[12,13]. The initial transverse momentum distribution 
9A(kT,Q 2 , b) is assumed to have a Gaussian form with a 
width that includes both an intrinsic part in a nucleon 
and nuclear broadening. Details of this model and sys- 
tematic data comparisons can be found in Ref. [9]. 

As demonstrated in recent studies, a direct conse- 
quence of parton energy loss is the medium modification 
of FF's [14,15] which can be well approximated by [16] 
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provided that the actual energy loss is about 1.6 times 
of the input value. Here z' c — Pt/(ptc ~~ A.E C ), 
z' g = (AL/X)pt/ AE c are the rescaled momentum frac- 
tions and AE C is the total parton energy loss during 
(AL/X) number of scatterings. The FF's in free space 
(z c , Q 2 ) are given by the BBK parameterization [17] . 
We assume a 1-dimensional expanding medium with a 
gluon density p g (r, r) that is proportional to the trans- 
verse profile of participant nucleons. The average number 
of scatterings along the parton propagating path is then 
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where z c = pt/ptc, V = He, <r(ab — * cd) are parton scat- where AL(b, r, 4>) is the distance a jet, produced at r, has 
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to travel along n at an azimuthal angle 4> relative to the 
reaction plane in a collision with impact-parameter b. 
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FIG. 1. Nuclear modification factors for hadron spectra in 
Au + Au collisions as compared to data from STAR [3] and 
PHENIX [2]. See text for a detailed explanation. 

According to recent theoretical studies [8,15,16] the to- 
tal parton energy loss in a finite and expanding medium 
can be written as a path integral, 
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p g (r,b,r + nT), (4) 



where po is the averaged initial gluon density at To in 
a central collision and (dE/dL)id is the average par- 
ton energy loss over a distance Ra in a 1-d expand- 
ing medium with an initial uniform gluon density pq. 
The corresponding energy loss in a static medium with 
a uniform gluon density po over a distance Ra is [16] 
dE /dL= (RA/2T )(dE/dL)id- In the high-energy limit, 
the parton energy loss has a logarithmic energy depen- 
dence [18]. However, for a parton with finite initial en- 
ergy the energy loss has a stronger energy dependence 
because of restricted phase space for bremsstrahlung [19]. 
Detailed balance between induced gluon emission and ab- 
sorption will further increase the energy dependence [20] 
for a parton with finite initial energy. In this paper we 
will use an effective quark energy loss 

(^>w = eo(E/n - 1.6) x -7(7.5 + E/p ), (5) 

from the numerical results in Ref. [20] . The detailed bal- 
ance reduces the effective parton energy loss and at the 
same time increases the energy dependence. The thresh- 
old is a consequence of gluon absorption that competes 



with bremsstrahlung and effectively shuts off energy loss 
for lower energy partons. 

Shown in Fig. 1 are the calculated nuclear modifica- 
tion factors Rab{pt) = da AB /(N^ )dcrb p for hadron 
spectra (\y\ < 0.5) in Au + Au collisions at y/s = 200 
GeV, as compared to experimental data. Here, (iV^n ) = 
J d 2 bd 2 rt A (r)tB(\b — r\). To fit the observed 7r° sup- 
pression in the most central collisions, we have used 
(solid lines) po — 1.5 GeV, eo = 1.07 GeV/fm and 
Ao = 1/(0730) = 0.3 fm with the new HIJING param- 
eterization [13] of parton shadowing. The hatched area 
(also in other figures in this paper) indicates a variation of 
6q = ±0.3 GeV/fm. The hatched boxes around Rab = 1 
represent experimental errors in overall normalization. 
Alternatively, one has to set po = 1.3 GeV and eo = 1.09 
when EKS parameterization [12] of parton shadowing is 
used (dot-dashed lines). Without parton energy loss, the 
spectra is slightly enhanced at px = 2 — 5 GeV/c due to 
nuclear kx broadening even with parton shadowing. 

The flat pt dependence of the 7r° suppression is a con- 
sequence of the strong energy dependence of the par- 
ton energy loss, which is also observed by other re- 
cent studies [21]. The slight rise of Rab at px < 4 
GeV/c in the calculation is due to the detailed balance 
effect in the effective parton energy loss. In this region, 
one expects the fragmentation picture to gradually lose 
its validity and other non-perturbative effects [22] be- 
come important that will give an enhanced (K + p) /% 
ratio in central Au + Au collisions. To include this 
effect, we add a soft component to kaon and baryon 
FF's that is proportional to the pion FF with a weight 
~ (NbinQ), r))/[l + exp(2pr c — 15)]. The functional form 
and parameters are adjusted so that [K +p)/ir sa 2 at 
Pt ~ 3 GeV/c in the most central Au + Au collisions 
and approaches its p + p value at pt > 5 GeV/c. The 
resultant suppression for total charged hadrons and the 
centrality dependence agree well with the STAR data. 
One can relate ft. ± and ir° suppression via the (K + p) /n 
ratio: R'X^ = R^ A [1 + {K +p)/n} AA /[l + {K +p)/tt] pp . It 
is clear from the data that (K + p) /tt becomes the same 
for Au + Au and p + p collisions at pr > 5 GeV/c. To 
demonstrate the sensitivity to the parameterized parton 
energy loss in the intermediate pr region, we also show 
R AA in 0-5% centrality (dashed line) for po = 2.0 GeV 
and eo = 2.04 GeV/fm without the soft component. 

Since jets produced in the central core of the dense 
medium are suppressed due to parton energy loss, only 
those jets that are produced near the surface emerge from 
the medium. The observed high-p^ hadron multiplic- 
ity should be proportional to the number of surviving 
jets in the outer layer of the overlapped volume which in 
turn is approximately proportional to the total number 
of participant nucleons. The nuclear modification factor 
normalized by (iVbi na ry) should decrease with centrality. 
This agrees well with the observed centrality dependence. 
This surface emission picture also gives a natural geomet- 
rical limit of the azimuthal anisotropy [23] . 
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In non-central collisions, the average path length of 
parton propagation will vary with the azimuthal angle 
relative to the reaction plane. This leads to an azimuthal 
dependence of the total parton energy loss and therefore 
azimuthal asymmetry of high-p^ hadron spectra [7,8]. 
Such asymmetry is another consequence of parton energy 
loss and yet it is not sensitive to the nuclear fcr broad- 
ening and parton shadowing. Shown in Fig. 2 is i>2(pt) 
(second Fourier coefficient of the azimuthal angle distri- 
bution) of charged hadrons generated from parton energy 
loss (dot-dashed) as compared to preliminary STAR data 
[24] using the 4-particle cumulant moments method [25] 
which is supposed to reduce non-geometrical effects such 
as inherent two-particle correlations from di-jet produc- 
tion [26] . The energy loss extracted from high-pT hadron 
spectra suppression can also account for the observed az- 
imuthal anisotropy at large pr- If the remaining v-i at in- 
termediate pt is made up by kaons and baryons from the 
soft component, one find that they must have w 2 ~ 0.23 
(0.11) for 20-50% (0-10%) collisions. The total v 2 (p T ) is 
shown by the solid lines. 
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FIG. 2. Azimuthal anisotropy in Au+Au collisions as com- 
pared to the STAR [24] 4-particle cumulant result. 

In the same LO pQCD parton model, one can also 
calculate di-hadron spectra, 

E ^ l^T* =^Y d 2 bd 2 rdx a dx b d 2 k aT d 2 k bT 

dz c dz d t A (r)t B (\h - r\)g A (k aT , r)g A (k bT , |b - r|) 
fa/A(x a , Q 2 , r)f b/B (xb, Q 2 , |b - r|) 



D h/C (z c , Q 2 , AE c )D h/d {z d , Q 2 , AE d ) 
s da 



2nz^z2 dt 



■(ab -> cd)S (p a +Pb~Pc- Pd), 
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for two back-to-back hadrons from independent fragmen- 
tation of the back-to-back jets. Set pn = P^' 8 , we define 
a hadron-triggered FF as the back-to-back correlation 
with respect to the triggered hadron: 



D^{z t ,^Pt S )=Pt 



trig 



d<7 A 1 A 2 /d 2 P T g dp T d(j) 
da h A \/d^ ' 



similarly to the direct-photon triggered FF [10] in 7- 
jet events. Here, zt = pt/Pt 1S and integration over 
|j/i,2| < Ay is implied. In a simple parton model, the 
two jets should be exactly back-to-back. The initial k^ 
distribution in our model will give rise to a Gaussian- 
like angular distribution. In addition, we also take into 
account the intra-jet distribution using a Gaussian form 
with a width of (k T ) = 0.8 GeV/c. 

Shown in Fig. 3 are the calculated back-to-back cor- 
relations for charged hadrons in Au + Au collisions as 
compared to the STAR data [4]. The same energy 
loss that is used to calculate single hadron suppres- 
sion and azimuthal anisotropy can also describe well the 
observed away-side hadron suppression and its central- 
ity dependence. In the data, a background B(pt)[1 + 
2v 2 (pt) cos(2A(j>)] from uncorrelated hadrons and az- 
imuthal anisotropy has been subtracted. The value of 
v%{j>t) is measured independently while B(pr) is deter- 
mined by fitting the observed correlation in the region 
0.75 < \cj)\ < 2.24 rad [4]. 
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FIG. 3. Back-to-back correlations for charged hadrons with 
p T rig > Pt > 2 GeV/c, p^ ig = 4-6 GeV/c and \y\ < 0.7 in 
Au + Au (lower curves) and p + p (upper curves) collisions as 
compared to the STAR [4] data. 



Since most of jets produced in the central core of the 
colliding volume are suppressed, the triggered high-p^ 
hadrons mainly come from jets produced near the sur- 
face traveling away from the dense core. These jets are 
not quenched, as observed in STAR data [4] . Most away- 
side jets, however, will be suppressed as they go through 
the dense core, except those that are propagating in di- 
rections tangent to the surface. This leads to a much 
smaller azimuthal anisotropy of the away-side suppres- 
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sion. On the average, both the magnitude of the away- 
side suppression and the centrality dependence should be 
similar to the single hadron suppression, as seen in the 
data. Contrary to another study [27], we find that the 
Pt broadening associated with energy loss has no sig- 
nificant effect on the observed back-to-back correlations, 
since those jets that have large final-state broadening also 
have large energy loss and thus are suppressed. 

Integrating over <^>, one obtains a hadron-triggered FF, 
D h ^{z Tl p^) = H l2 d^D h ^{z T AJT^)- Shown in 
Fig. 4 are the suppression factors of the hadron-triggered 
FF's for different values of p^ ls in central Au + Au 
collisions as compared to a STAR data point that is 
obtained by integrating the observed correlation over 
7r/2 < |A0| < 7r. The dashed lines illustrate the small 
suppression of back-to-back correlations due to the ini- 
tial nuclear broadening in d + A collisions. The strong 
QCD scale dependence onp^ lg of FF's is mostly canceled 
in the suppression factor. The approximately universal 
shape reflects the weak px dependence of the hadron 
spectra suppression factor in Fig. 1, due to a unique en- 
ergy dependence of parton energy loss. 
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FIG. 4. The suppression factor for hadron-triggered frag- 
mentation functions in central (0-5%) Au + Au (d+Au) colli- 
sions as compared to the STAR data [4]. 

In summary, we have studied simultaneously the sup- 
pression of hadron spectra and back-to-back correlations, 
and high-pT azimuthal anisotropy in high-energy heavy- 
ion collisions within a single LO pQCD parton model 
incorporating current theoretical understanding of par- 
ton energy loss. Experimental data of Au + Au col- 
lisions from RHIC can be quantitatively described by 
jet quenching in an expanding medium. With HIJING 
(EKS) parton shadowing, the extracted average energy 
loss for a 10 GeV quark in the expanding medium is 
{dE/dL) ld re 0.85(0.99) ± 0.24 GeV/fm, which is equiv- 
alent to dE /dL pa 13.8(16.1) ± 3.9 GeV/fm in a static 
and uniform medium over a distance Ra = 6.5 fm. This 
value is about a factor of 2 larger than a previous esti- 
mate [16] because of the variation of gluon density along 
the propagation path and the more precise RHIC data 
considered here . 
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